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Abstract. Superconductivity, superfluidity, condensation, cluster formation, 
etc. are phenomena that might occur in many-particle systems. These are due 
to residual interactions between the particles. To explain these phenomena con- 
sistently in a microscopic approach, at some point, one needs to solve few-body 
equations that are modified because of the Pauli principle and the interactions 
of the many particles around. 

1 Introduction 

Correlations in many-particle systems are responsible for a number of interest- 
ing and exciting properties of the system. In an infinite, equilibrated system 
the presentation of these properties (or phases of matter) are usually organized 
in a phase diagram. This is a plot of the dominant state of matter separated by 
"critical lines" depending on the two basic variational parameters of quantum 
statistics relevant for a grand canonical ensemble, i.e. the temperature T and 
the chemical potential /i. Instead of the chemical potential that is difficult 
to access experimentally one might by use of proper thermodynamic relations 
chose the particle density n{fi,T). Note, however, that the relation between fj. 
and n{^,T) is not simple and depends on the specific system considered. 

The phase diagram becomes particularly rich when the chemical potential fi 
is in the order of MeV to GeV. This is at rather large particle densities, where 
nuclei are close together. This is the domain where Quantum Chromodynamics 
(QCD) dominates the dynamics of the system. Since nucleons are made of 
quarks it is obvious that the relevant degrees of freedom to describe the system 
might change from nucleons (nuclear matter) to quarks (quark matter) at some 
point. In addition gluons might appear as additional degrees of freedom that are 
otherwise absent at lower densities. Because of the free (color) charge this phase 
is called a (quark gluon) plasma. Indeed, lattice QCD and model calculations 
suggest a variety of phenomena in the QCD dominated phase of matter. These 
are superfluidity, critical points, first and second order phase transitions, color 
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Figure 1. Schematic view of the phase diagram of nuclear matter. The phase diagram 
is empirical accessible by heavy ion collisions, excited nuclei, observation of neutron 
stars and the early universe as indicated in the diagram. New plans at GSI aim at 
exploring the color superconducting phase as well. 



super conductivity, plasma phase, among others. A sketch of the phase diagram 
is given in Fig. ^ To explain all these different states of matter one needs a 
treatment that goes beyond the simple quasiparticle picture. 

To tackle the complicated many-particle system the Green function ap- 
proach is a good starting point 0]. Equations for Green functions are hierar- 
chally coupled and hence an approximate solution of the problem is unavoid- 
able. If we are interested in two-body correlations only, the three-body Green 
functions are approximated. Quite a few phenomena are already accessible at 
this level. However, there are good reasons to go beyond two-particle correla- 
tions, e.g.: 

• Particle production even in a dense environment such as deuteron for- 
mation in a heavy ion reaction, need a third particle to conserve energy- 
momentum [Q. 

• To study the properties of a-particles j^] or determine the critical tem- 
perature of a possible a-particle condensate ^ needs an in-medium 
four-body equation. 

• Recent results in the Hubbard model indicate, that three-particle con- 
tributions may lead to a different (lower) critical temperature compared 
to the simple Thouless criterion m. Question of this type have not been 
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addressed for nuclear matter. 

• The chiral phase transition is often discussed along with a confinement- 
deconfinement transition based on investigating mesons (quark-antiquark 
states). Does this transition happen for nucleons (three-quark states) at 
the same density /temperature? 

Some of these issues related to the nuclear matter phase are addressed in the 
present paper. The ones related to quark matter are given by Mattiellio et 
al[|. 

2 Theory 

We use Dyson equations to tackle the many-particle problem, see e.g. Ref. §. 
This enables us to decouple the hierarchy of Green functions. The Dyson equa- 
tion approach used here is based on two ingredients: i) all particles of a cluster 
are taken at equal time ii) the ensemble averaging for a cluster is done for an 
uncorrelated medium. The resulting decoupled Green functions may be eco- 
nomically written as resolvents in the 71-body space, where n = 2, 3, 4, ... is the 
number of particles in the considered cluster. 

The solution of the one-particle problem in Hartree-Fock approximation 
leads to the following quasi-particle energy 

^1 = ^ + E ^2(12,T2)/2 ^ -4ff + ^""(0). (1) 
2mi ^ 2mf 

The last equation introduces the effective mass that is a valid concept for the 
rather low densities considered here and = 11 — S^^{0). The Fermi function 
fi = f{ei) for the z-th particle is given by 

^(^^) = e.(.-.) + l - (2) 
The resolvent Go for n noninteracting quasiparticles is 

n 

G^{z) = {z~H^)-^N = Ro{z)N, H^^Y.''- (3) 

where Go, i?o, and N are formally matrices in n particle space. The Matsubara 
frequency z\ has been analytically continued into the complex plane, za 2: 0]. 
The Pauli-blocking for n-particles is 

iV = /l/2.../„±/l/2.../„, /-I-/ (4) 

where the upper sign is for Fermi-type and the lower for Bose type clusters. 
The full resolvent G{z) is given by 

G{z)^{z-H^-V)-'N, y^Y. ^^V^- (5) 

pairs a 
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Note that ^ V. For the two-body case as weh as for a two-body subsystem 
erabedded in the n-body cluster the standard definition of the t matrix leads 
to the Feynman-Galitskii equation for finite temperature and densities Q], 

T^i^) = + V,"N^Roiz)T^{z). (6) 

Introducing the Alt Grassberger Sandhas (AGS) transition operator Uap{z) 
the effective inhomogeneous in-medium AGS equation reads 

U^p{z) = (1 - 5^p)R^\z) + N]T^{z)Ro{z)U-,p{z). (7) 

The homogeneous in-medium AGS equation uses the form factors defined by 

\Fp) ^Y.~^0,N]V,^\i;B,) (8) 

7 

to calculate the bound state ipBa 

\Fo.) = J2^c,pN^Ti{B3)Ro{B3)\Fp). (9) 

Finally, the four-body bound state is described by 

= E SarU;^{B4)Ro{B4)N]T^{B,)R,{Bi)\T:;), (10) 

where a C u, 7 C r and ct, r denote the four-body partitions. The two-body 
input is given in (^) and the three-body input by (^. Note that, although 
we have managed to rewrite the above equations in a way close to the ones 
for the isolated case, they contain all the relevant in-medium corrections in 
a systematic way, i.e. correct Pauli-blocking and self energy corrections. The 
numerical solution requires some mild approximations that are however well 
understood in the context of the isolated few-body problem. 

3 Results 

An experiment to explore the equation of state of nuclear matter is heavy ion 
collisions at various energies. Here we focus on intermediate to low scattering 
energies and compare results to an experiment ^■^^Xe-t-^^^Sn at 50 MeV/A by 
the INDRA collaboration |0. A microscopic approach to tackle the heavy 
ion collision is given by the Boltzmann equation for different particle distribu- 
tions and solved via a Boltzmann Uchling Uhlenbcck (BUU) simulation ||]. 
The reaction rates appearing in the collision integrals are a priori medium 
dependent. We use the in-medium AGS equations (^ that reproduce the ex- 
perimental data in the limit of an isolated three-body system. For details on the 
specific interaction model see Ref. . We investigate the influence of medium 
dependent rates in the BUU simulation of the heavy ion collision as compared 
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Figure 2. BUU simulation of deuteron Figure 3. Ratio of proton to deuteron 

formation in the collision of ^^^Xe+^^^Sn numbers as a function of cm. energy. The 

at 50 MeV/ A. Use of in- medium rates lead experimental data are from the INDRA 

to a 20% enhancement. collaboration. 



to use of isolated rates. Figure g shows that the net effect (gain-minus-loss) 
of deuteron production becomes larger for the use of in-medium rates (solid) 
compared to using the isolated rates (dashed). The change is significant, how- 
ever, a comparison with experimental data is difficult since deuterons may also 
be evaporating from larger clusters that has not been taken into account in 
the present calculation so far. The ratio of protons to deuterons may be better 
suited for a comparison to experiments that is shown in Figure |^. The use of 
in-mcdium rates (solid) lead to a shape closer to the experimental data (dots) 
than the use of isolated rates (dashed). Within hnear response theory for infi- 
nite nuclear matter the use of in-medium rates leads to faster time scales for 
the deuteron life time and the chemical relaxation time as has been shown in 
detail in Refs. jl^, This faster time scales should have consequences for the 
freeze out of fragments. 
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Figure 4. Due to the Mott effect clusters 
exist only above the respective curves. Be- 
low no bound states exist because of Pauli 
blocking. 



Figure 5. Part of the phase diagram with 
lines of critical temperatures. Solid line 
four-body AGS type calculation. Others 
earlier calculation Ml 
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As does the nucleon, see eq. (|l|), the cluster changes its mass. The binding 
energy changes as well until the Pauli blocking is too strong for bound states to 
exist (Mott effect). This effect depends on the momentum of the cluster. This 
important effect for modeling of heavy ion collisions is shown in Figure ^. 

In Figure |^ part of the phase diagram of nuclear matter is shown. The 
condition for the onset of superfluidity for a-particles is B(Tc, /ic, P = 0) = 4/ic, 
where B is the binding energy. The critical temperature found by solving the 
homogeneous AGS equation for ^ < confirms the onset of a condensation 
even at higher values (solid line) as given earlier (dashed, from jij) which was 
based on a variational calculation using the 2+2 component of the a particle. 
For yti > the condition E — Afi for the phase transition can also be fulfilled. 
However, the significance for a possible quartetting needs further investigation. 

To conclude, strongly correlated Fermi systems such as nuclear/quark mat- 
ter provide an excellent field for few-body techniques. 
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